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A WIREBASKET PRECONDITIONER FOR THE 
MORTAR BOUNDARY ELEMENT METHOD 


THOMAS FtiHRER AND NORBERT HEUER 

Abstract. We present and analyze a preconditioner of the additive Schwarz type for the 
mortar boundary element method. As a basic splitting, on each subdomain we separate 
the degrees of freedom related to its boundary from the inner degrees of freedom. The 
corresponding wirebasket-type space decomposition is stable up to logarithmic terms. For 
the blocks that correspond to the inner degrees of freedom standard preconditioners for the 
hypersingular integral operator on open boundaries can be used. For the boundary and 
interface parts as well as the Lagrangian multiplier space, simple diagonal preconditioners 
are optimal. Our technique applies to quasi-uniform and non-uniform meshes of shape- 
regular elements. Numerical experiments on triangular and quadrilateral meshes confirm 
theoretical bounds for condition and MINRES iteration numbers. 


1. Introduction 

In recent years, different variants of the non-conforming boundary element method (BEM) 
have been developed. The underlying boundary integral equation is of the hrst kind with hy¬ 
persingular operator. Non-conformity refers to the presence of discontinuous basis functions. 
(Note that, in the case of integral equations of the second kind or hrst kind equations with 
weakly-singular operator, conforming basis functions can be discontinuous.) The hrst paper 
on non-conforming BEM considers a Lagrangian multiplier to deal with the homogeneous 
boundary condition on open surfaces jS]. This technique was extended in [HI to domain 
decomposition approximations, and is usually referred to as mortar method. In this paper 
we study preconditioners for the mortar BEM presented in [H]. These are the hrst results 
on preconditioning techniques for linear systems stemming from non-conforming boundary 
elements. 

Our preconditioner is based on a decomposition of the approximation space and choosing 
locally equivalent bilinear forms. It therefore hts the additive Schwarz framework. There is a 
large amount of literature on the additive Schwarz method, mainly aiming at hnite element 
systems, see, e.g., dSl EHl [21] for overviews. For additive Schwarz techniques applied to 
boundary elements dealing with hypersingular operators see, e.g., [mEiiES], cf. also |22) 
for an overview. In particular, m considers a wirebasket-oriented splitting. Graded meshes 
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on curves, locally refined and anisotropic meshes (on surfaces) have been analyzed, respec¬ 
tively, in [3 [ini [16]. Other variants, also for hypersingular operators, consider overlapping 
decompositions and multiplicative applications, see, e.g., [25l[n]. 

In this paper we extend the additive Schwarz technique to mortar boundary elements. 
In this case, due to the presence of a Lagrangian multiplier, system matrices have a saddle 
point structure. This structure can be handled by using standard arguments aiming at the 
minimum residual method (MINRES). More precisely, the spectrum of the system matrix 
(with or without preconditioner) is being controlled by the spectrum of the main block and 
the singular values of the off-diagonal block (arising due to the presence of the Lagrangian 
multiplier), see |2^. A second complication due to the non-conformity of the method is that 
the bilinear form a(-, •) representing the hypersingular operator is replaced by the weakly sin¬ 
gular operator acting on surface differential operators (surface curl). There are no standard 
preconditioners for this bilinear form. Our strategy is to split the subspace of discontinuous 
basis functions from the rest of the approximation space X^. The remainder X^ forms 
a subspace of the energy space of the hypersingular operator. It turns out that the bilinear 
form a(-, •) reduces to the standard one of the hypersingular operator when restricted to 
Xj^. In this way, standard preconditioners like the ones mentioned previously can be applied 
to this block (actually, there are individual blocks associated with each subdomain). Now, 
the other subspace X^ contains all the basis functions associated to interface or boundary 
nodes. In domain decomposition terms, it is a wirebasket space and is related with the skele¬ 
ton (or wirebasket) of a coarse mesh which is formed by the subdomains of the underlying 
decomposition. In our case, basis functions associated to the boundary of a subdomain T, 
can be decoupled from the other elements of X^, they form a subspace Xj i. It turns out 
that the bilinear form a(-, •) restricted to Xj i is spectrally equivalent to a diagonal matrix 
(a mass matrix related to the boundary of T*). In this way, simple diagonal matrices can 
be used (for the preconditioner) to reduce the problem of preconditioning the bilinear form 
a(-, •) : Xh X Xh —)■ M to the standard one of hypersingular operators on each subdomain. 
In our numerical examples we will use multilevel diagonal scaling from j3| for these parts. 

A priori error analysis for domain-oriented non-conforming boundary elements yields 
quasi-optimal error estimates which are perturbed by (poly-) logarithmic terms depending 
on the mesh size, see mm- These perturbations appear due to the non-existence of a well- 
defined trace operator in the energy space of hypersingular operators. It is unknown whether 
estimates of these perturbations are sharp. Naturally, such logarithmic perturbations also 
appear in the analysis of additive Schwarz preconditioners, at least when considering non¬ 
overlapping decompositions. Note that, in our method, we subtract a wirebasket space and 
this amounts to trace operations at the boundaries of subdomains. Also in the case of finite 
elements, such splitting operations cause logarithmic perturbations, see, e.g., [5]. Even¬ 
tually, our main result considers combinations of simple diagonal and multilevel diagonal 
preconditioners and proves that they are optimal up to poly-logarithmic terms. In some 
cases, logarithmic perturbations of condition number bounds can be optimized by multiply¬ 
ing terms of the preconditioner by different logarithmic weights, see again, e.g., jl]. In this 
paper we consider three different weightings {Cases 1,2,3) where Cases 2,3 are optimized to 
show a bound 0(|log(h)|^) for the condition number of the preconditioned system. Here, h 
denotes the minimum of the diameters of all elements. In contrast, for Case 1 (which does 
not use logarithmic weights for the parts dealing with the bilinear form a(-, •)) the theoretical 
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bound is 0(|log(/i)|^), worse than the bounds for Cases 2,3. Our numerical experiments, on 
the other hand, indicate that Case 1 is superior to Cases 2,3. This suggests that some of 
the theoretical bounds used for the analysis are not sharp, at least in the particular situation 
of our numerical examples. 

An outline of the remainder of this paper is as follows. In the next section we present the 
model problem, recall the dehnition of some Sobolev norms, and present a mortar discretiza¬ 
tion for the model problem. In Section [3] we recall some results on the MINRES method, 
present our subspace decompositions and corresponding preconditioners, and state the main 
result (Theorem [9]). Proofs are given in Section 13.31 Some numerical experiments are re¬ 
ported in Section 01 They all conhrm our theoretical estimates from Theorem |9l though 
exhibit smaller logarithmic perturbations than predicted. In particular, we also study the 
case of locally rehned meshes driven by adaptivity, where preconditioners behave as expected. 
Let us note that we do not know of any a posteriori error analysis for mortar boundary el¬ 
ements. The only known results concerning non-conforming BEM consider two-level (or 
h — h/2) estimators applied to the Nitsche coupling |0j, not the mortar coupling. 

Notation. We abbreviate estimates of the form A < C ■ B with some constant C* > 0 by 
A < B. In particular, we use this notation if C is independent of the mesh size and the 
number of elements. Analogously, we use A> B for A > C ■ B. If both A < B and A> B 
hold true, we use the notation A B. Moreover, |x| denotes the Euclidean norm for a point 
X e M^. 


2. Mortar boundary elements 

In this section we briefly recall some results on the mortar boundary element method. 


2.1. Model problem and functional analytic setting. Let T C x {0} denote a 
plane open surface with polygonal boundary ST. For simplicity we refer to T as a domain 
in 

We recall some dehnitions of Sobolev spaces. Let S' C be a bounded subset and dehne 
for 0 < s < 1 the seminorm 


u 


2 

H^iS) 



u{x) - u{y)\^ 
\x — |/|2(«+1) 


dx dy. 


Then, H^{S) is equipped with the norm 

Illl2 ||||2 Ii2 


and H^{S) is dehned as the completion of C^{S) with respect to the norm 




H^S) + dist(a:, dSy^ 


dx, 


where dist(a:, 5S') := inf \x — y\. The dual spaces of H^{S), resp. iL®(S'), are denoted by 

^ y£dS 

resp. H~^{S). Additionally, (•, ■)s denotes the L^(S) scalar product, which is 

continuously extended to the duality pairing on x H^{S), resp. x H^{S). 
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Let n G denote a normal vector on L, e.g., n = (0,0,1)^. Define the hypersingular 
integral operator (formally) by 

Wm{x) := - A /■„(;,) 

oUa, Jr dUy \y - x\ 

It is well known that this operator extends to a continuous mapping between and 

/f-V2(r). 

Our model problem reads as follows: Given / G -L^(r) we seek for a solution u G 
such that 

(1) {Wu,v)r = {f,v)r Wv E 

The usual conforming boundary element method consists in replacing by a finite¬ 
dimensional subspace Xh C and seeking for a solution Uh G Xh such that 

{Wuh , Vh)r = if , Vh)r Vh/i G X^. 

In this paper, we study preconditioners for a non-conforming scheme that is based on a 
decomposition of the surface T. In the next section we introduce the corresponding subspace 
decomposition. The non-conforming method based on this decomposition is called mortar 
method and is presented in Section 12.31 

2.2. Subspace decomposition and meshes. Let Ti,..., T^v denote a decomposition into 
non-intersecting (open) polygonal subdomains giving rise to the coarse mesh 

N 

r:={ri,...,rjv} with r = ljr,-. 

i=i 

Each subdomain T* is equipped with (a sequence of) regular and quasi-uniform meshes 
The minimum and maximum diameters of elements of the meshes 71 are denoted by and 
hi, respectively. As in [U] we assume without loss of generality that h, < 1 and set 

h := min h, and h := max A. 

Let JCi denote the set of nodes of We will also need the set of interior nodes /C° and the 
set of nodes on the boundary of T,, X} := /C, \ /C°. For a node Zj G Xi we denote by the 

(bi)linear basis function which satisfies rl'j\zk) = Sjk for all Zk G /Cj. Introducing the space 
of piecewise (bi)linear functions 

Xh,i := {n G G°(ri) : V = ^ with aj G M} 

Zj 

we define the product spaces 

N N 

Xh := HXh, C := 

i=l i=\ 

We denote the respective degrees of freedom by Ki := jj^Xi = dim(A'/j j) and K := — 

dim(A/i). Note that by definition of Xi, elements of Xh do not necessarily satisfy the homo¬ 
geneous boundary condition on c^T nor continuity across interfaces dTi fl dTj. 

4 











r 





^ 4 







Ti 

Yi2 









^34 












lis 






r, 










Ti 








Ti2 

Yi3 










^2 


T 23 








Figure 1. Subspace decompositions with non-conforming meshes. 

For the mortar BEM, we denote the interface of two neighboring subdomains Fj 7 ^ F^ by 
'jij := int(Fj n F^) where “int” refers to the (relative) interior. Also, let diam(S') denote the 
diameter of a set S' C We will need the following assumption. 

Assumption 1. Each non-empty interface 'jij = consists of an entire 

edge of Vi or F^. //diam(5F fl c^Fi) > 0, then 5F fl (9Fj is a union of edges o/Fj. 

Figure [T] shows examples of two different subspace decompositions with non-conforming 
meshes. Given the skeleton 

N 

7:=U^r„ 

i=l 

we infer from Assumption [ 1 ] that 7 is covered by a set of non-intersecting interface edges and 
boundary edges 

L 

r := { 7 i,..., 7 l} with 7 = ^ 7 ^. 

i=i 

For each interface edge 7 ^, let £iag G N} resp. i^or G N} denote the indices 

of the subdomains F^j^^ resp. such that 

7 ^ = 77 ag,^mor aud 7 ^ is an edge of 

In particular, we set £iag := i if 7 ^ is a boundary edge, i.e., 7 ^ C c^F, n 5F so that we can 
handle the homogeneous boundary conditions and the interface conditions simultaneously. 
On each 7 ^ we introduce a mesh such that the following assumption is satisfied. 

Assumption 2. The mesh 0 / 7 ^ is a strict coarsening In particular, any 

element of Ti covers at least two elements 

Moreover, the diameters of the elements in T£ and are comparable, i.e., there exists 

a eonstant C > 0 such that 

C~^ diam(f) < diam(T) < diam(f) for all t E T£ and T G with T G t. 

The constant C > 0 is independent of 1. 
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We also define the discrete spaces on edges, 

Yh,£ ■= {"0 e ■ "01* is constant for all t G r^}, i = 1,..., L, 


and the (global) space for the Lagrangian mnltiplier 

L 

Yh-.= l[Yh,£. 

i=i 

Based on the previous definitions of decompositions, meshes, and subspaces, we next intro¬ 
duce the mortar boundary element method. 


2.3. Mortar BEM. Let v = (ui,... jVn) with sufficiently smooth component functions vj 
defined on Lj. We define the piecewise differential operator curlji/ by 


N 


cuAhv := y^(curlr,nt)° with curlr-n* = {dyVi{x, y), -d^Vi{x, y), 0) 


1=1 


and (•)° being the extension by 0 onto L. 

We need the single layer integral operator 

V^(x) : = 


\x-y\ 


dy, 


;:-V2 


_ -j /o 

which extends to a continuous operator, mapping (L) to (L). Here 

:= {(ni,n2,n3) G {HY\T)f : Ug = 0} 

.—'- 1/2 

and (L) is its dual space. Furthermore, we define the jumps [n] across interface edges 

7£ by 


|«]|„ := _ 

With the definitions of the bilinear forms 


if 7£ C dT, 
else. 


N 


a[u,v) := (Hcurl^M, curlj^n) 7 - ;= ^^(Hcurl/fM, curlrin)ri. 


2=1 


b{u,'ip) := ([n] , '0)r := , ' 0 ) 7 £ 

£=1 

for all u,v E and 'ifj G for some e > 0 and the right-hand-side functional 

N 

Y{v) := 5 ^(/, Vi)r,, 


2=1 


we can state the mortar BEM: Find {uh, 4>h) G x W such that 

a{uh,Vh)+ b{vh,(l)h) = F{vh), 

b{uh,i>h) = 0 


( 2 ) 


for all {vhiiJh) G Xh X W. This formulation admits a unique solution. 
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Theorem 1 ([9l Theorem 2.1]). Let Assumptions [IHl hold true. There exists a unique 
solution {uh, (ph) £ 'xYh of ([2]) . Assume that the exact solution u of satisfies u G 
/fi/^+''(r) for some r G (0,1/2], Then, there holds 

\\u Ufi\\H^/'^{T) ~ |lc>S(Zl)| h ||m|| 

Remark 2. The work [9] deals with homogeneous boundary conditions on dV. However, 
the analysis can be generalized to the present situation, see [8] for the case of BEM with 
Lagrangian multipliers. In particular, [9l Theorem 2.1] holds true if we do not impose ho¬ 
mogeneous boundary conditions in X^. 

We note that the bilinear form a{-, ■) is not elliptic due to the fact that piecewise constant 
functions c = (ci,..., Cn) G Xh with q G M are in the kernel of curl//(-). Therefore, for our 
analysis we will use a simple stabilization of a(-, •) which is similar to the one that is often 
used for hypersingular integral equations on closed surfaces. 

Lemma 3. Let 0 7 ^ a G M denote an arbitrary but fixed constant and let G W denote the 
characteristic function on 7 ^, i.e., = dtk o-nd define 

L 

(3) a(M, n) := a(M, n) + for all u,v ^ X^. 

i=i 

Then, the variational equation ([2]) is equivalent to: Find {uh, 4>h) ^ Xh x W such that 

a{uh,Vh)+ b{vh,(ph) = F{vh) 

biuhy'iph) = 0 

for all {vh,'ilJh) eXh-xYh. 

Moreover, foru,v G Xh, there holds 

(5) ]|n]|^i/ 2 (^) < |log(h)|a(n,n) and a{u,v) < \log{h)\^\\u\\Hi/^r)MHi/ 2 (^r)- 
The involved constants do not depend on h. 

Proof. To see the equivalence, we note that b{uh,^e) = 0 from the second equation of ([2]) 
resp. (jl]), since G Yh. Hence, the additional stabilization terms in the definition of a(-, •) 
always vanish. Note that the solutions {uh, fJh) of (E]) and (jl]) are in fact identical. 

The upper bound in ([ID follows from the continuity 

a{u,v) < |log(h)| 2 ]|M]|j:^i/ 2 (r)l|w]|Hi/ 2 (r) 
of a(-, •) (see m Lemma 3.9]) and the continuity 

biu^fj) < \\og{h)\Y^\\u\\Hi/2^r)M\L^{A 
of fe(-, •) (see P Lemma 3.14]), since 

L L 

< |log(h)|]|n]|^i/2(r)lkl|j^i/2(r)5^l|l|| 

e=i i=i 

To derive the lower bound in (jSD we note that the analysis from [ 8 ] and P yields 

N 

(6) l^l^i/2(r) •= ~ fo’^ all ^ ^ 

i=l 
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Moreover, we apply the following result from [9l Proposition 3.5], which comes from a discrete 
Poincare-Friedrichs inequality for fractional-order Sobolev spaces proved in [10]: There exists 
a constant (7 > 0 such that for all e G ( 0 , 1 / 2 ] and any v G with u|ar = 0 there 

holds 


(7) 


C iu]|i2(r) < e >|^i/2+e(r) + diam(7£) 

7^ is interior edge 
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Let P D P denote an extension of P with (9P C P. Moreover, let T denote a subdomain 
decomposition of P with T C T such that the shape regularities of T and T are equivalent. 
In particular, T can be chosen such that each boundary edge C c?Pn5Pj is an interior edge 
in T. Thus, ([^ holds true if we replace P, resp. T, with P, resp. T. For each v G 
we set v\-j- = V and ul^- 7 - := 0. Then, v G and llw]|j:^i/ 2 +e(;y) = ]|u] 7 i/ 2 +e(- 7 -). We 

infer that 


( 8 ) 




^lr/i/2+qf) 


L 

^diam(7£)“^“2^ 

i=i 


L 

= £"VlHi/2+7r) + 

e=i 



2 


L 

e=i 


with some constant C > 0 depending on a and the diameters of 7 ^ but not on e. Finally, 
choosing e = |log(h)|“^, the inverse estimate 


4^|u7i/2+e(7-) < \v\m/2(r) 

together with ([HD and ([HD shows ellipticity of a(-, •). 


□ 


2.4. Discretizations. For real-valued vectors we use bold symbols, e.g., x. Each vector 
X G is uniquely associated to a function v G Xh in the following way. Let 


bnd') 

I'/l ) ■ ■ ■ ) 'Ikv ''1 ’ • • • ’ 'IK21 • • • ’ '/l ) ■ ■ ■ ) 'IKm J 


db ^(2) 


( 2 ) 






denote the basis of Xh- For simplicity we use the notation pi,... ,? 7 i^ for the basis. Then, 
X G (with K = Tdj) corresponds to 


N Ki K 

-EE ■> = E ^jVj- 

i=l j=l j=l 

We define the Galerkin matrix A G of a(-, •) as 


Ajk := a{r]k, r]j) for j,k = l,..., K. 


Let {xf'^} denote the basis of with xf'^ltk = ^jk for tk G T£. Analogously as before, 
we write Xi, ■ ■ ■, Xm with M := ■— for the corresponding basis of W. 


Then, each if} ^Yh can be written as 


M 

■0 = YjXj foi' some y G 
i=i 

We define the matrix B E by 


Bjk ;= b{r]k, Xj) j = I, ■ ■ ■, M, k = 1,..., K. 


Denoting the right-hand side vector by f G with := F{j]k)-, the formnlation (jTj) is 
eqnivalent to the matrix-vector equation; Find (x, y)^ G such that 


C 



A 

B 0 




3. Preconditioning 

In this section we analyze different wirebasket preconditioners for the mortar BEM con¬ 
sidered in Section [2J First, we recall results on the MINRES method. 


3.1. Minimal residual method. Throughout we consider the preconditioned minimal 
residual method (MINRES) with inner products (x, y)p ;= y^Px induced by block-diagonal 
preconditioners of the form 



where the blocks Pa G Pb G are symmetric and positive definite. (Here 

and in the following, empty spaces represent null matrices of appropriate dimensions.) The 
preconditioned system then reads 


p-^c 


Furthermore, define the matrix 


IPb^B ) 


fA B^\ __ /P^^/^AP^P^/^B^Pb^/^ 

Vb J ■“ l^Pg'/^BP^'/' 


p-l/2cp-V2_ 


We note that there holds spec(C) = spec(P ^C) for the respective spectra. Let 
(10) Amin < Amax and El < • ■ ■ < Em 


denote, respectively, the extremal eigenvalues of A and the nonzero singular values of B. Of 
course, they are all positive. We also define the condition number 

k{C) := max{|A|; A G spec(C)}/min{|A|; A G spec(C)}. 

The following is a well-established result, see, e.g., |27| . 
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Proposition 4. Denote by := P ^(f — the residual of the k-th preconditioned 

MINRES iteration with inner product {■, •)p. Then there holds 


^ Up < 2 

r(o)||2, - ^ 


^(C)-A 

k(C) + 1 J 


so that the number of preconditioned MINRES iterations, required to reduce the initial resid¬ 
ual to a certain percentage, is bounded by 0(k(C)). 


Bounds for the spectrum of C can be specified in terms of the eigenvalues of A and singular 
values of B. 


Proposition 5 ( |191 Lemma 2.1]). There holds 

spec(p-‘C) C|i(A„,„ - \/AL„ + 4Ef„, 1(A„„ - y'Af.„ + 4E?)] 

u [Amin! 2(^max + \/Amax + 


with Amin? Amax? being the numbers from (ITU]) . 


3.2. Preconditioner and main results. Our preconditioning technique is based on an 
initial decomposition of into wirebasket components related with the coarse mesh T and 
the remainder. Then, individual preconditioners are applied to the three spaces of wirebasket 
and interior components and the Lagrangian multiplier. 


3.2.1. Wirebasket splitting. For the initial decomposition of Xh, we dehne for each Vi G Xh^i 
the unique representation 

(11) Vi = Vi^i + Vifi with Vifl G Xifl := {w G X^^i : wlan = 0} and Vi^i := Vi - Vifl, 


and the preconditioning forms dj : X^ x Xf, —?■ M (j = 1, 2, 3) dehned by 


N 


N 


(12a) 

di{u, v) : 


9r, 

■}L^dTi) 

+ ^{WiUi^o , Wi,o)ri, 



N 



1 ^ 

(12b) 

d2{u,v) : 

= '^{UipldFi , 
2=1 

ar, 


1 1, 2;..M9 y , Vi,o) 

|log(h)|2 



N 



N 

(12c) 

di,{u,v) : 

= 5^|log(h)|(ni,i|ar, 

! ar, 

:)i,2(ari) + ^ ]{WiUifl , nj,o)ri 


2=1 2 = 1 


i 5 


Here, Wi is the hypersingular integral operator associated with the subdomain Fp Note that, 
with Vi being the simple-layer integral operator associated with Fj, we have 

(l/iCurlr,ni,o! curlr,ni,o)ri = {WiVifl, Vi^o)ri ^ lbi,o|||i/2(r,)- 

Dehnition (IT2l) provides, up to logarithmic terms, stable splittings. 
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Lemma 6. For all v G Xh there holds 

\^og{h)\-^di{v,v) < a{v,v) < \\og{h)\‘^di{v,v), 
\^og{h)\-^d 2 {v,v) < a{v,v) < \\og{h)\'^d 2 {v,v), 
\^og{h)\-^d 3 {v,v) < a{v,v) < \\og{h)\d 3 {v,v). 

A proof of Lemma [H] is given in Section 13.31 


3.2.2. Preconditioner for A. We now consider a preconditioner for the matrix A that corre¬ 
sponds to the bilinear form a(-, •) on Xh x Xh- Having performed the initial decomposition 
of Xh into wirebasket and interior components, Lemma [6] and the structure of the bilinear 
forms dj defined by flT^ show that it suffices to provide preconditioners for the L‘^{dri) 
terms and the terms involving the hypersingular integral operator. We use, respectively, a 
simple diagonal preconditioner and an arbitrary preconditioner for the hypersingular integral 
operator in the conforming case, see, e.g., ^mmm- 

In the following, let 'Pwi denote such a preconditioner for the hypersingular integral op¬ 
erator Wi with constants Amax such that 

(13) < {WiVifl , < A^^x^Pw^^x 


for all Vifi G q with Uj o = Furthermore, let PaPi denote a preconditioner 

with 


(14) 


/^mLy^ParW < |IWil9rJli2(arp < hmLy^P^rA 


for all Vi^i G Xi^i with Vi^i = • Define the preconditioner Pd) g for the 

f-th subdomain by 


p('0 



Pw/i/ 

\log{h)\-^Pw. 


/|log(h)|Par, 


V 


Wi 


if di(-, ■) is used, 
if d 2 (-, ■) is used. 


if dsl', ■) is used. 


and the overall preconditioner Pa for the matrix A, corresponding to the bilinear form a(-, 
on Xh, by 


/p{i) 


Pa := 


\ 


•) 


For the last two definitions we have assumed an appropriate order of the degrees of freedom in 
Xh^i- The logarithmic terms in the definition of P*^*^ stem from the logarithmic perturbations 
in the definition (I12p of dj(-, ■). In the remainder of this work we will refer to 
‘'Case j” if dj{-, •) is used in the definition of Pa (j = 1, 2, 3). 

Our main result concerning the preconditioning of A is as follows. 
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Theorem 7. Set 


Kin ■= ..., and A^ax := max{A^L) /^mL) • • •, 

Then, there holds for all x G 

( |log(/i)|“^AminX^PAX < x^Ax < |log(/i) ^Amaxx'^PAX for Case 1, 

I |log(/i)|“2AminX^PAX < x^Ax < |log(/i) ^Amaxx'^PAX for Case 2, 

[ |log(/i)|“^AminX^PAX < x^Ax < | log(/i) | AmaxX^P ax for Case 3. 

Therefore, the condition number of A is hounded by 

«(A) < \log(h)fK 

'^min 

with (3 = 5 in Case 1 and f3 = 4 in Cases 2,3. 

Proof. The proof follows directly from Lemma Eland Assumptions fll3D . (1141) on the precon¬ 
ditioners. □ 

3.2.3. Final preconditioner for the full matrix C. In order to dehne the preconditioner P (jO]) 
for the full matrix C we assume that we have a matrix Pb G such that there exist 

numbers amm, o'max > 0 with 

(15) ^rainy^Pey < ||V'lli 2 (.^) ^ ^max y^Pny 

for all fj = ymXm G Yh- Below, we will select Pb to be diagonal with or without 

logarithmic scaling. 

To provide bounds for the spectrum of C by means of Proposition [5] it remains to bound 
the singular values Si,..., S^ of the matrix B. 

Lemma 8. Let 0 < Si < • • ■ < S^ denote the nonzero singular values of the matrix B and 
let Amin, Amax defined as in Theorem\f^ Then, 

{ Amin<^min|log(h)| ^ Sj^ ^ ^ Amax^max f^^ CaSe 1 , 

AminCrmin|log(h)|"^ < S? < S^ < AmaxC^max for CaSe 2, 

AminCrmin|log(h)|"2 < S? < S^ < AmaxC^max|log(h) |“^ for CaSe 3. 

A proof of Lemma [ 8 ] will be given in Section 13.31 

Now, let M G denote the ^^( 7 ) mass matrix, i.e., 

^jk ■= {Xj : Xk)'y for j,/c = 1,...,M. 

Obviously, M is diagonal and 

M 

= y^My for all V' = X] XmXm e Yh- 

m=l 

The main result of our paper is the next theorem. Its proof is immediate by combining the 
previously established estimates, namely Theorem [7] and Lemma [HI together with the general 
results provided by Propositions ID and 0 

Theorem 9. Let Amin, Amax defined as in Theorem\^ and let (Tmin, o'max > 0 6e the numbers 
from (lT5l) . Then the spectrum of the preconditioned matrix has a superset like spec(C) C 
[—a, —b] U [c, d] with numbers a,b,c,d > 0 that satisfy the following estimates. 
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• Case 1; //Pb = |log(/i)| then cimin = o-max = |log(4)| and 

Ainin/max{Ama^, Aj^/f^}|log(/i)|-=^ <b<a< , 

Amm|log(/i)r^ <c<d< max{Amax, A^^fx}|log(4)P- 

• Case 2; //Pb = |log(/i)|~^M, then cimin = o'max = |log(4)| and 

Amin/max{Amax, Ai/f^}|log(/i)|"^ <b<a< X]llJlog{h)\^/^, 

Amm|log(4)r^ <c<d< max{Amax, A^^^} |log(4) 

• Case 3; //Pb = M, then cTmin = o'max = 1 and 

Ai„in/max{Amax, Aj„/f^}|log(/l)|“^ <b<a< (Amax)^^^|log(/^)^^/^ 

Amin|log(/i)r^ <c<d< max{Amax, A^/f^}|log(/i)|. 

Therefore, the condition number of C is bounded by 

«(C) < |log(/i)|^max{Aj^/f^,A^^^}/Amin 

with (d = 5 in Case 1 and fd = A in Cases 2,3. Furthermore, the number of preconditioned 
MINRES iterations, required to reduce the relative residual to a certain threshold, is bounded 
like the condition number in the respective case. 

3.3. Proofs and technical details. For the proof of Lemma E] we need a trace inequality 
and an inverse estimate, which are given in the following two lemmas. 

Lemma 10 ([H Lemma 4.3]). Let R (Z Mf be a bounded Lipschitz domain. There exists a 
constant C > 0 such that for all e G (0,1/2) holds 

ll^^l|L2(aij) < for all v G //^/^+^(/?). 

Lemma 11 ([121 Lemma 4]). For a function Vi G Xi^h with splitting fITT]) . Vi = Vifi + Vi^i, 
there holds 

||wi,o|liJi/2(r,) < |log(4i)lllwi||Hi/2(rq, i = l,...,N. 

The proof of |121 Lemma 4] uses [51 Lemma 4.5]. An alternative proof of Lemma ITT] which 
utilizes multilevel norms is given in [TTl Theorem 3.6]. 

Proof of Lemma\^ We start with a proof of the upper bound. Let v = G 

with n® := Vi^o and := Vi^i. Application of the triangle inequality, boundedness (jS]) of 
the bilinear form a(-, •), and equivalence (V ■, •)r — ]| • ]|?_ 1 / 2 ^ ^ together with the estimate 

]| ■ ]|— 1/2 , < ]| ■ ]|— 1/2 for fractional-order Sobolev spaces, leads to 

Ht (r) (T) 

a{v,v) < F + |log(h)n]n(^)]|^i/ 2 (^) 

N N 

< ]|curlr,ni,ol|~-i/2^^^^ + |log(h)P ]|ni,i]|^i/2(r,) 
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Then, {Vi- , •)ri ^ || • ||~-i/ 2 , , and (l^jCurlriMi,o, curlrini,o)ri = {WtUi^ , nj,o)ri for all G 

(rd 

Xh show 


N N 

a(y,'y) < '^{WiVi^o, ni,o)r, + \\og{h)\'^'^\\vi^i\\ly2(^r^y 

i=l i=l 


Let Tj denote a closed extension of the subdomain Tj and let 71 denote an extension of 
the mesh Ti such that the shape-regularities of the meshes 71 and % are equivalent. For 
Zj G /C° we set := and for Zj G X} we define as the (bi-)hnear function with 
vf\zk) = Sjk for all nodes Zk of the mesh 71- Hence, an arbitrary function 

Vi = 'Yhzj&Ki ^ define its extension Vi as 

-Zj 


By the properties of the and i7^/^-norms, we have 

II l |2 < ||~ Il 2 

\\Vi,l\\H^/2{Y.) ^ ||'«^i,l||^i/2(rp- 

Set ujk := supp(^*^). We note that there exists a constant Cco\ > 0 that depends only on the 
shape-regularity of the mesh 71 such that 

lixfchfc 

ZfcG/Cj 


With 


11^7 ||2_ 

W'lk lljyl/2(^^) 


- diam(a;fc) ~ ||h?|9rjli2(^^n9rd = 
and the locality of the L^-norms we further deduce 


^11^1/2(0;^) ~ i 11 Vk n^r^) 




iWl'^idri)- 


Thus, altogether we have 


N N 

a{v,v)<^{W,v„,,v,,)r. + \\og(h)\^Yl ||'^i,i|9rillL2(ari)) 

i=l i=l 

which proves the upper bounds. 

For the lower bounds, we use Lemma fTOl with i? = Fj and e = |log(h)|~^ (for h small 
enough). Together with an inverse inequality this gives 


Using the norm equivalence 

{WiVifl , 


lWi|9rJli2(arp ^ |fog(4)|||u||^i/2(r,). 

II ■ 111/1/2 (r ) — (Wi' ) ■)ri and Lemma fTTl shows that 

'r'i,o)ri — ||'f^i,o||^i/2(-p.) ^ |fog(if)l ll'^i|li/i/2(rp- 
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Combining the previous relations and summing over i = 1,..., N proves 

N N 

di{v,v) = 'yi,o)r, + 5^||wi.il9rJli2(9rp < \^og{h)\^\\v\\\^/ 2 ^ry 


2 = 1 
N 


2=1 


N 


d2{v,v) = ^|log(h)| ^(PhiUi,o, 'yi,o)r, + Iki.i|arj|i2(arp < \^og{h)\\\v\\l,^2^^^, 


2 = 1 
N 


2=1 


N 


d 3 {v,v) = Ui,o)ri + |log(h)| ^ WviAdTiWh^aVi) < \^og{li)\^\\v\Wy 2 (jy 


2 = 1 


2 = 1 


Hence, by applying the ellipticity of a(-, •) from TheoremlH this shows the lower bounds. □ 

Proof of Lemma O Note that the nonzero singular values of B are given by the square roots 
of the eigenvalues of the matrix BB^, since B has full (row) rank. Furthermore, we note 
that the smallest and largest singular values are given, respectively, by the minimum and 
maximum of the term 

K M 


— with V yki]k and ^ x^Xm- 

y€i.x ||y||p.,l[x|[p, ^ 


max 


We start with the upper bound. By the Cauchy-Schwarz and triangle inequalities we have 

N 

< 2||V'||i2(.^)^ ||ui||i2(arp < l|V'|li2(^)|log(4)|^^dj(u, u) 


2=1 


with Ml = M 2 = 0 and M 3 = —1. This together with dj{v,v) < AmaxY^PAY and < 

tTmaxX^PRX from flT^ proves the upper bound. 

For the lower bound, we use the proof of Lemma [ 6 ] to see that 

AminY^PAY ^ dj{v,v) < |log(h) ||u|| 

with rrii = m 3 = —2 and m 2 = —1. This leads to the estimate 


max -—— > (Amin)^'^^|log(h)| max - 

yGR^\{0} ||yIIpaI|x||pb ^ 6 A;,\{o} ||n||m/ 2 (r)II^H pb 


6(u,'0) 


By using ||' 0 |li 2 (..y) > CTminX^PBX from (fT5]) and the discrete inf-sup condition 

6(u,'0) 


sup 


11 II 


> / 3 ||V'||l 2 (^) for all fj &Yh 


from jH Lemma 3.12] (with constant /9 > 0 independent of h) we conclude the lower bound. 

□ 


4. Numerical Examples 

In this section we present numerical examples in which we compare the different behaviors 
of the preconditioned systems induced by the preconditioning forms di{-, ■), d 2 {-, ■), d 3 {-, ■). 
Note that block Pa of the preconditioner P is determined by the preconditioning forms 
dj{-, •). For the second block Pb we choose, up to a possible logarithmic term, the (diagonal) 

15 






mass matrix M for the Lagrangian multiplier space. We distinguish the following cases (with 
corresponding numbers (Tmin, o’max according to Theorem |9l cf. flT^ b 


Case la) 

PB = iiog(/i)riM, 

<^min ~ <^max ~ |lo§(A)|5 

Case lb) 

Pb = M, 

^min ^max I 5 

Case 2) 

Pb = |log(/i)|-iM, 

<^min ~ <^max ~ |log(/i)|, 

Case 3) 

Pb = M, 

^min ^max !• 


Note that Case la, Case 2, Case 3 correspond to the bounds obtained in Theorem[9l whereas, 
at least theoretically, we would expect worse bounds for Case 1 b. Moreover, we compare the 
results to a simple diagonal preconditioner with 


(Pa)^^ ^jj^jk Slid (Pb)^^ ^jj^jky 


where 6jk denotes the Kronecker delta symbol. In the figures and tables below, we refer to 
this preconditioner as diag. 

Throughout, we use the MINRES algorithm, see Section l3Tl to solve the discrete system. 
We stop the algorithm if the relative residual in the fc-th step satisfies 


,(fc)| 


~( 0 )| 


< 10 


-6 


4.1. Diagonal preconditioner and mnltilevel diagonal preconditioner. For the 

wirebasket component we use a simple diagonal preconditioner. Indeed, it is straightforward 
to prove that 

(16) Ik*,i|arjli2(gr,) ~ ^ ^ 

zjeK] 


where yjvf ^ ^ diam(a;j) is the diameter of the node patch Uj = 

supp(h]*^) ^j- Ihe example from Section we define the diagonal preconditioner 
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djk- 


According to (IT^ it is optimal, that is, the numbers from (114)) behave like 
(17) /i^)^ ~ pW, ~ 1, i = 


For the example from Section IT3] we test as preconditioner (for the wirebasket components) 
the diagonal of the matrix, i.e., we set 

(ParJifc := {V^cur\HVk^ , cuvlHvf^)ri5jk. 

Since (l/jCurlji/p)*^, cur\Hrij^'^)ri — diam(ci;j) ~ fhe constants from (114)) satisfy (ITT)) in 

this case as well. 

It remains to select preconditioners for the matrix blocks that belong to the interior 
unknowns, i.e., the ones corresponding to the nodes /C°, i = 1,... ,N. As indicated by 
()T3)) . it is enough to take for each subdomain a standard preconditioner that works for the 
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hypersingular operator. In the following we use as Piy. a multilevel diagonal preconditioner, 

i.e., 


p-i 


£i=0 


D7^T 




More precisely, we consider 7i = 71,Li as the hnest level of a sequence of meshes 7i/ {i = 
0,..., Li). Then, D^. is the diagonal part of the Galerkin matrix of (PcurljL(-), curljL(-))ri 
with respect to the nodal basis of Xj,o on level ii and T^, is the matrix representation of the 
embedding operator which embeds elements of the space Xj,o on a coarse level ii to functions 
on the hne level Lj. For the examples from Section we replace the entries of by 

|(^q|i/ 2 /i 2 _ Here, cuj* is the support of the basis functions of level ii associated with node j. 

It is known, see, e.g., [T], that these preconditioners are optimal on triangular meshes, i.e., 
the constants from (IT^ satisfy 

(18) — 1 
\ 'mm 'max 


with mesh size independent constants. Such multilevel preconditioners can be extended 
to locally rehned meshes with assumptions on rehnement zones, see, e.g., jT], or by use of 
special rehnement strategies like Newest Vertex Bisection, cf. [7]. The basic idea is that 
smoothing with the diagonal elements is done with respect to the degrees of freedom, where 
the associated basis functions have changed. 

We remark that the cited results for the multilevel diagonal preconditioners are stated for 
triangular meshes only. However, for uniform rehnements, the same techniques can be used 
to prove optimality on quadrilateral meshes. Finally, note that (fT7|) and (1T8|) imply that the 
numbers Amin, Amax from Theorem [7] satisfy 


A • ~ A ~ 

^'min — /'mfl.x — 


1 . 


According to Theorem |9] we then expect bounds k(C) = 0(|log4|^) with /3 < 5 in Case 
la and < 4 in Cases 2,3. A theoretical bound for the condition number in Case lb 
would results in an exponent (3 > 5 (and is not given here). But our numerical results show 
that this preconditioner is as competitive as in Case la, and better than in Cases 2 and 3. 
Our explanation is that some of the technical bounds used in proofs are not sharp, see the 
discussion in the introduction. 


4.2. Problem on Z-shaped domain with triangular meshes. We consider the vari¬ 
ational formulation (jl]) with / = 1 and stabilization parameter a = 0.1 on the Z-shaped 
domain F from Figure[2l In this case we consider only one subdomain, i.e., T = {F}, V = 1. 
The stabilization parameter is chosen such that the lower order stabilization terms in the 
dehnition ([3]) are not the dominating parts in the condition numbers (for large h). 

For the dehnition of the Lagrangian multiplier space W, we combine two adjacent bound¬ 
ary edges to one element of the mesh r^. Mesh rehnement is driven by Newest Vertex Bisec¬ 
tion, see, e.g., [TB]. In particular, we note that each triangle T is divided into 4 son elements 
Ti,... ,Ti, with \Tj \ = |T|/4. Moreover, this rehnement rule preserves shape-regularity, i.e. 


sup 

TeTi 


diam(T)^ 

|r| 


S sup 

TgTo 


diam(T)^ 

in 
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Figure 2. Initial triangulation of Z-shaped domain for the example from Section 1121 
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Figure 3. Condition numbers of the preconditioned systems and number of 
iterations in the MINRES algorithm for the example of Sect ion 112] and uniform 
rehnements. 


which also holds for adaptive mesh rehnements. For details we refer the interested reader 
to [15] and references therein. We remark that the initial triangulation does not satisfy 
Assumption [2] since, for instance, the boundary edge (—1, 0) x {0} x {0} of To contains only 
one boundary element. We use a uniform rehnement in the hrst step, which ensures that 
Assumption [2] holds true. 

Figure m shows the condition numbers (left) as well as the numbers of iterations (right) 
needed to reduce the relative residual in the MINRES method by 10“® in the case of uniform 
rehnements. Additionally, the condition numbers are listed in Table H 
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step 

dof 

h 

k 

no prec. 

la 

lb 

2 

3 

diag. 

1 

49 

1/2 

1/2 

25.19 

6.09 

6.52 

8.03 

6.44 

6.60 

2 

153 

1/4 

1/4 

26.33 

7.97 

7.51 

9.56 

8.31 

10.77 

3 

529 

1/8 

1/8 

46.68 

9.70 

8.66 

27.66 

13.50 

21.95 

4 

1953 

1/16 

1/16 

94.46 

11.22 

9.76 

59.18 

21.41 

45.98 

5 

7489 

1/32 

1/32 

192.36 

12.60 

10.77 

104.18 

30.11 

97.20 

6 

29313 

1/64 

1/64 

392.15 

13.88 

11.75 

162.97 

39.24 

206.00 

7 

115969 

1/128 

1/128 

799.48 

15.12 

12.77 

235.40 

48.56 

436.41 


Table 1. Condition nnmbers of the preconditioned systems for the example 
of Section 14.21 and nniform rehnements. 
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Figure 4. Logarithmic behavior of condition nnmbers for example of Sec¬ 
tion I4.2I and nniform rehnements. 


In the following let ns refer to Pj as the preconditioner of “ Case j” {j G {la, 16, 2, 3}). The 
nnmerical results indicate that the preconditioners Pia and Pib sue better than the others, 
and that P 3 is better than P 2 . In contrast. Theorem [9] predicts better bounds for P 2 and 
F 3 . Nevertheless, all the results conhrm the theoretical estimates. Indeed, Figure H] indicates 
that k(C) is bounded by 0(|log(h)p) even in Case 2. 

In the next example we consider adaptive mesh rehnements, where we use a simple ZZ- 
type estimator, see, e.g., [ 6 ], to mark elements for rehnement and additionally rehne all 
elements that share a boundary edge. We note that this estimator is not analyzed in [ 6 ] 
for the present (non-conforming) situation, but is heuristically used to obtain adaptively 
rehned meshes. Condition numbers of the preconditioned systems and numbers of iterations 
needed in the MINRES algorithm are plotted in Figure [5l Moreover, the condition numbers 
are listed in Table O We observe similar results as in the case of uniform rehnements. In 
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Figure 5. Condition numbers of the preconditioned systems and number 
of iterations in the MINRES algorithm for the example of Section I4.2I and 
adaptive rehnements. 


step 

dof 

h 

h 

no prec. 

la 

lb 

2 

3 

diag. 

1 

49 

5.00e-01 

1/2 

25.19 

6.09 

6.52 

8.03 

6.44 

6.60 

2 

98 

5.00e-01 

1/4 

22.51 

6.09 

5.73 

8.17 

6.40 

6.88 

3 

202 

5.00e-01 

1/8 

23.05 

7.09 

6.46 

26.00 

12.58 

8.49 

4 

444 

5.00e-01 

1/16 

23.74 

8.37 

7.91 

57.98 

20.98 

11.35 

5 

939 

5.00e-01 

1/32 

44.60 

9.40 

8.69 

95.06 

27.59 

14.42 

6 

1961 

3.54e-01 

1/64 

75.78 

10.31 

9.42 

158.09 

38.06 

18.36 

7 

4038 

2.50e-01 

1/128 

107.28 

11.21 

10.71 

245.60 

50.67 

24.21 

8 

8289 

1.77e-01 

1/256 

168.65 

12.13 

11.62 

345.67 

62.42 

30.52 

9 

16939 

1.25e-01 

1/512 

240.84 

13.20 

12.46 

470.05 

75.43 

43.72 

10 

34516 

1.25e-01 

1/1024 

408.58 

14.28 

13.15 

612.52 

88.46 

57.62 

11 

70278 

8.84e-02 

1/2048 

755.93 

17.34 

15.54 

882.78 

115.87 

79.32 


Table 2. Condition numbers of the preconditioned systems for the example 
of Section I4.2I and adaptive rehnements. 


particular, our theoretical results are conhrmed also for adaptively rehned meshes. Again, 
the results for the weakest of the domain decomposition preconditioners, P^., indicate that 
^ 0(|log(h)|^) also in this case, cf. Figure |6l 

4.3. Problem with four subdomains and quadrilateral meshes. We consider the 
variational formulation (|3]) with / = 1 and stabilization parameter a = 0.1 on the quadratic 
domain F := (0,2)^ x {0} with a decomposition into four subdomains Fi = (0,1)^ x {0}, 
r 2 = (1,2) X (0,1) X {0}, Fa = (0,1) x (1,2) x {0}, and r 4 = (1,2)^ x {0} sketched in 
Figure [71 For the intersections Fi n r 2 and Fi fl Fa, we dehne Fi to be the Lagrangian side 
and for the intersections r 4 fl r 2 and r 4 Cl Fa, we dehne r 4 to be the Lagrangian side. We 
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Figure 6. Logarithmic behavior of condition numbers for example of Sec¬ 
tion |1]2] and adaptive refinements. 


define the Lagrangian elements that come from Fi and r 4 as the union of two adjacent edges 
that lie in dTi. For the Lagrangian elements that come from r 2 and Fs we take the union of 
three adjacent boundary edges. 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

X 


Figure 7. Subspace decomposition of F = (0,2)^ x {0} and their initial 
meshes for the example from Section I4.31 
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Figure 8. Condition numbers of the preconditioned systems and number 
of iterations in the MINRES algorithm for the example of Section 14.31 with 
different rehnenient levels of subdomain meshes. 


step 

dof 

hi 

Iki 

hs 

^4 

no prec. 

la 

lb 

2 

3 

diag. 

1 

172 

1/4 

1/6 

1/6 

1/4 

34.77 

13.97 

12.94 

17.72 

13.96 

17.43 

2 

360 

1/4 

1/12 

1/6 

1/8 

63.62 

17.67 

15.95 

43.41 

21.77 

27.83 

3 

948 

1/4 

1/24 

1/12 

1/8 

194.24 

21.10 

18.70 

85.29 

31.43 

50.64 

4 

2804 

1/8 

1/48 

1/12 

1/8 

446.57 

25.27 

22.39 

156.85 

46.10 

97.47 

5 

10532 

1/8 

1/96 

1/24 

1/16 

1545.48 

32.08 

28.76 

279.10 

69.07 

204.58 

6 

39492 

1/16 

1/192 

1/24 

1/32 

3661.95 

40.65 

37.04 

476.84 

101.74 

436.13 

7 

45316 

1/32 

1/192 

1/48 

1/64 

1948.27 

47.65 

43.99 

572.20 

122.29 

477.40 


Table 3. Condition numbers of the preconditioned systems for the example 
of Section 14.31 with different rehnement levels of subdomain meshes. 


We consider uniform rehnements where each element of Tj is divided into four elements. 
For the experiment we rehne each of the subdomain meshes separately, which leads to differ¬ 
ent mesh sizes hi, h 2 i lhi Ih- Note that for our problem conhguration there holds hj = hy 
The results are given in Figure [ 8 ] and Table |3l As in Section [4.21 we observe that the pre¬ 
conditioners Plai Pib corresponding to the preconditioning form di(-, •) behave best in terms 
of condition numbers and numbers of iterations. The preconditioners P 2 and P 3 stemming, 
respectively, from d 2 {-, ■) and d 3 (-, •) show a stronger dependence on the mesh size. Neverthe¬ 
less, theoretical bounds are conhrmed also for this example. In particular. Figure [3 suggests 
that k(C) ;< 0(|log(h)p) for all domain decomposition preconditioners. 

Let us also remark that the condition number of the un-preconditioned system gets smaller 
from Step 6 to Step 7, see Table [3l The condition number k(C) is bounded (up to logarithmic 
terms) by Amax/Amin and this term depends on the ratio jh. Since h gets smaller and h 
stays constant from Step 6 to Step 7 (as we rehne all subdomains except r 2 ), this explains 
the observation. 
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Figure 9. Logarithmic behavior of condition numbers for example of Sec¬ 
tion |1]3] with different refinement levels of subdomain meshes. 
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